Dispersion Relations for Bernstein Waves in a Relativistic Pair Plasma 
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A fully relativistic treatment of Bernstein waves in an electron-positron pair plasma has remained 
too formidable a task owing to the very complex nature of the problem. In this article, we perform 
contour integration of the dielectric response function and numerically compute the dispersion curves 
for a uniform, magnetized, relativistic electron-positron pair plasma. The behavior of the dispersion 
solution for several cases with different plasma temperatures is highlighted. In particular, we find two 
wave modes that exist only for large wavelengths and frequencies similar to the cyclotron frequency 
in a moderately relativistic pair plasma. The results presented here have important implications for 
the study of those objects where a hot magnetized electron-positron plasma plays a fundamental 
role in generating the observed radiation. 
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I. INTRODUCTION 



Relativistic electron-positron pair plasmas are found 
in many astrophysical objects such as neutron star mag- 
netospheres fl], relativistic jets and accretion disks asso- 
ciated with black holes in the centers of active galactic 
nuclei 0, 01 ■ To study the properties of the pair plasma 
in these objects it is imperative to employ a fully rela- 
tivistic approach. In the study of magnetized pair plas- 
mas one finds that a number of oscillation modes exist 
(see [1, §)• The focus of this study is to investigate the 
behavior of Bernstein modes in a uniform, magnetized, 
relativistic e'^e^ pair plasma. Bernstein waves are elec- 
trostatic undulations that are always localized near the 
electron cyclotron harmonics in an electron-ion plasma 
. These waves can propagate undamped only very close 
to the plane perpendicular to the static magnetic field. 
Such waves are of great interest since, in an electron-ion 
plasma, they strongly interact with the electrons and are 
excellent candidates for plasma heating and driving cur- 
rents as compared with the electromagnetic (O)rdinary 
and the e(X)traordinary modes 

The non-relativistic treatment of Bernstein waves in 
an electron-ion plasma is well understood [1, d, [l^, [HI ■ 
It is the fully relativistic case that is marred by difficul- 
ties such that a closed form analytic solution is hard to 
formulate or maybe even impossible. Many workers have 
expounded on the fully relativistic treatment of electron 
Bernstein waves [12], the ultrarelativistic case [l3|, and 
have successfully obtained approximate dispersion rela- 
tions [ll,[il[i6j. The Bernstein mode in a weakly rela- 
tivistic pair plasma has been investigated by the authors 
of [l3| where they have found closed curve dispersion 
relations that are remarkably distinct from the classical 



case. All of these studies either discuss the fully relativis- 
tic case to the point where no closed form analytic solu- 
tion is found, or simplify the analysis by either treating 
the limiting case only or employ various approximations 
that may not yield an entirely correct result. 

In this paper, we present a fully relativistic treatment 
of Bernstein waves in a pair plasma and provide disper- 
sion curves that highlight the transition from the weakly 
to strongly relativistic regime. The rest of the article is 
organized as the following. We derive the key equations 
for the dielectric response tensor and describe our nu- 
merical approach to the problem in the following section 
(§IIIl- In Section Hm we present the main results of this 
study along with a discussion on how the solution behaves 
as a function of the plasma temperature, and plasma fre- 
quency. In the final section (§ lIVp we highlight some of 
the important points of the study. 



II. RELATIVISTIC DISPERSION RELATION 

The evolution of the distribution function, /(r, p, t), of 
plasma particles in phase space is governed by the Vlasov 
equation, which in the momentum representation is given 
as 



dt 



v-Vr/s+g.(E + vxB)-Vp/, -0 (1) 
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where s indicates different species constituting the 
plasma. To investigate the behavior of small amplitude 
waves with oscillation periods much smaller than particle 
collision times, we make the following assumptions 

/,(r,p,i) = /o.(p) + /i.(r,p,t) (2) 
B = Bo + Bie^(''--"*' (3) 
E = Eie'(*'-'"~'^*) (4) 

where the subscripts and 1 indicate equilibrium and 
perturbed functions, respectively. Furthermore, we re- 



2 



strict the equilibrium distribution function to only de- 
pend on the momentum of the particles to account for 
the anisotropy introduced by the ambient static mag- 
netic field. As a result, we write the Vlasov equation in 
its linearized form 

jj,s = -g.(Ei + V X BOe^C^— *) • Vp/o. (5) 

The main idea here is to calculate the perturbation of 
the distribution function by integrating along the unper- 
turbed orbits from some time in the past, say to the 
present time t. Next, we use Ohm's law to write the 
current density induced by the perturbed distribution 

s s 

This enables us to write the effective dielectric permit- 
tivity tensor in terms of the conductivity tensor a Q 



dfojp) 

dpii 



X ■ — : ■ p sm 



sm TT^a 



•mk\ 



7 dfo{p) 
P± dp± 



p^smOdpdO (11) 



In the above, || and _L subscripts denote components par- 
allel and perpendicular to the equilibrium magnetic field, 
J-ya{0 is the Bessel function of non-integer order, J!ya{0 
is the derivative of the Bessel function with respect to ^, 
where ^ = ^^^5^, and a = ^ with Uc denoting the non- 
relativistic cyclotron frequency. Finally, one finds the 
dispersion relation, to — uj(k), by setting the dielectric 
response function to zero. 



e(w,k) =k-e(u;,k) - k^O 



(12) 



e{uj,k) = eo M - 



lUJseo 



(7) 



In the fully relativistic approximation, the energy and 
linear momentum of particles in the rest frame of the 
plasma are given as 



p — 7mv 



(8) 



where the Lorentz factor is given in terms of the momen- 
tum as 



7= 1 



p- 



(9) 



The complete details of the rest of the calculation can be 
found in various monographs and textbooks on plasma 
waves (see [l,[i3), and we only provide the salient points 
of the derivation in what follows. We adopt Bq = Bqz for 
the equilibrium magnetic field and restrict both the wave 
vector and the perturbed electric field to be k = A;_Li and 
El = Eix as dictated by the purely electrostatic mode 
where the perturbed magnetic field Bivanishes. After 
some mathematical manipulations we arrive at the rela- 
tivistic dielectric tensor 



e('^,k) = 








(10) 



where the different components have been summed over 
both species, e"*" and e~, of the pair plasma 



1 



sin irja 



7 dfoip) 2 
-1 j p sm 



p_L dp±_ 



P± dfoip) 



ujmeQ J LUc dp± [_ sin ir^a 



Xj^aiOJ-^aiO + ^ }p' Sin edpdd 



which in our case simply picks out the exx component. 
To keep the treatment fully relativistic we adopt the 
Maxwell-Boltzmann- Jiitner distribution function iH, UM , 



.Up) - (47rm^c^) 



K2{V) 



where 



(13) 



(14) 



is the ratio of the rest mass energy of the particles to that 
of their thermal energy, and K2 is the modified Bessel 
function of the second kind and of order two. Also, the 
equilibrium plasma distribution has been normalized to 
unity 



(15) 



l = no= / h{p)d-'p 



Taking the derivative of /o(p) with respect to p±_ and 
yields 



9p^i> 



(16) 



where <|b can be replaced by either || or _L components. 
At this point, we can carry out the integration over the 
polar angle and by defining /3 — 
the following 



we can write exx as 



2 2 



k\nv'c'K2{r}) Jq 



p^e-"^ 



(17) 



X 



■J.yaif3 sm0)J^^aiP sin0) — 1 > sin 6d6dp 



!q |^sm7r7a 

where the non-relativistic plasma frequency is defined as 

=2 



UJ. 



p 



»-oe 



(18) 
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Next we use the following Bessel function identity Q 

sm0Ja{bsme)J^a{bsine)de (19) 







2 sin 7ra ^ / 1 3 , , 

2F3 -, 1; -, 1 - a, 1 + a; -b^ 

ira \ Z 2 



to express the integral over the polar angle in terms of a 
hypergeometric function, and redefine all constants and 
variables to make e^x dimensionless 



P 

mc 



k± = 



UJQ 



(20) 



LUr, 



— (3 — [3 = k±p a 

UJQ 



U! 



Then, we can write exx as the following where the integral 
now is just over the dimensionless momentum p, 



kl 

2' ' 2' 



V 



(21) 



71!;, 1 + 7c2); — /3 ) dp — 1 



In writing this equation we have made use of the following 
integral identity 



p e ' 'dp — 

V 



(22) 



The integrand in Eq. (|2ip consists of a hypergeometric 
function which is singular for 1 — "fuj — —n for integer 
n = 0, 1, 2, . . . ,. This singular behavior, as we shall see, 
is associated to the phenomenon of cyclotron resonance 
where plasma particles are in resonance with the wave at 
the cyclotron harmonics. Also, this very resonance poses 
a real challenge for any numerical computation of the in- 
tegral and has to be dealt with using advanced numerical 
techniques. Since our interest lies in finding the oscilla- 
tion frequency w as a function of the wavenumber /cj^, 
it is clear from Eq. (|2ip that this operation is explicitly 
nonlinear. Therefore, one is left with an exercise of root 
finding for a given a). Alternatively, one could simplify 
the analysis by making some approximation. However, 
by adopting such methodology one risks losing the sub- 
tleties of the solution and may obtain something that is 
not entirely correct, as we show in the weakly relativistic 
case. We remain optimistic and decide to compute the 
dispersion relation using a brute force method, that is by 
simply integrating Eq. (PT|). 



A. Numerical Approach 

As the Lorentz factor is a function of momentum, the 
integrand remains singular over the domain of integra- 
tion. A workaround for avoiding the singular points on 
the real axis is by analytically continuing the momentum 



to the complex domain. We can easily shift the integra- 
tion contour below the real axis by writing p p ~ i5 
where 5 is reasonably small. Ideally, one would like to 
keep the contour on the real axis but go below the sin- 
gular points to avoid divergence while following the Lan- 
dau prescription [19;]. A similar result can be achieved by 
closing the contour of integration in the lower half of the 
complex plane as shown in Fig. [T] Here we are interested 



Im p 



Re p 



FIG. 1: Contour of integration in the complex p-plane used 
to avoid singular points of the integrand. Only the general 
case is shown here to guide the reader. 



in finding the principal value of the integral in Eq. (|2ip 
which is given by the following, 

/"OO 

p.V. / dp f{p) ^ITT ^0 - + /3 + h) (23) 

^ Residue 



where 



h = 
h = 



dp fip) « 



-iS 



dp f{p) 

00 — iS 


dp f{p) 



(24) 



and Rq is the residue from the poles on the Re(p) axis. 
In the above, for I2 we note that the integrand vanishes 
sufficiently rapidly for large p due to the exponential. 
When Lu and k± are kept real, one finds that Rq is also 
real which makes the first term on the right in Eq. (j23p 
purely imaginary. However, the second term has both 
real and imaginary components, and the negative sign is 
indicative of the clockwise sense of the contour. Since the 
left hand side is real then so must be the right, which sug- 
gests that the imaginary components cancel each other 
and yields the following result 



p.V./ #/(p) = -Re(/3+/4 



(25) 



Although it appears that with the given prescription 
one can have an arbitrarily large S, we find that the inte- 
gral does start to lose its accuracy as S approaches unity. 
Therefore, we set S = 0.1 for all numerical computations. 
To compute the integral numerically over the hypergeo- 
metric function we used Mathematica (V.6) for it is capa- 
ble of calculating generalized hypergeometric functions. 
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The poles of the integrand were dealt with by employ- 
ing a globally adaptive integration routine available in 
Mathematica. To speed up the integration over the sin- 
gular points we used the Double Exponential Quadrature 
singularity handler built into Mathematica's integration 
routine. 



III. PROPERTIES OF DISPERSION CURVES 

In the non-relativistic case of the electron Bernstein 
modes in an electron-ion plasma, one finds that there 
are no wave modes below the first harmonic. The dis- 
persion curves above the hybrid frequency are all bell 
shaped with local maxima corresponding to stationary 
modes. Furthermore, band gaps are present above the 
hybrid frequency between each dispersion curve. It is 
not at all surprising to say that the picture is remarkably 
different in the relativistic pair plasma scenario. We plot 
the dispersion curves for a relativistic pair plasma in Fig. 
[21 [31 [4] for different values of 77 — 1, 5, 20, respectively. We 
assume a plasma frequency of a)p = 3 for these plots. 

The dispersion relation for a moderately relativistic 
pair plasma, shown in Fig. [2l clearly has two wave 
modes. This is further accompanied by the existence 
of two stationary modes with vanishing group velocity 
[duj/dk^ — 0) at two distinct oscillation frequencies for 
fcj^ = 0. Also, above the higher stationary mode there are 
two wavenumber solutions for a given w. This behavior 
persists in the case of a mildly relativistic pair plasma, 
shown in Fig. [3l However, one sees some drastic changes 
in the shape of the dispersion curves as the particles lose 
their energy. We readily notice the appearance of the 
curve near the cyclotron fundamental frequency. This 
marks the onset of the cyclotron resonance where the 
plasma particles oscillate at the same frequency as the 
perturbing electrostatic wave. Although not very signif- 
icant at this point, the higher harmonic resonances also 
start to emerge. Furthermore, the dispersion relation 
now extends to higher wavenumbers and the turnover 
from the lower wave mode into the upper mode is not as 
sharp as it was in the previous case. We also notice a 
shift in the stationary modes, and we discuss this point 
in a later section. 

The authors of [T^ reported the dispersion relation 
for a weakly relativistic pair plasma. They found island 
shaped curves occurring mostly between the cyclotron 
harmonics. We find a similar solution for the weakly rel- 
ativistic case, shown in Fig. [H but with a few exceptions. 
Firstly, the dispersion curves are no longer closed but ex- 
tend to higher wavenumbers as they approach cyclotron 
harmonics. This behavior is very similar to what we ob- 
serve in the case of a non-relativistic electron-ion plasma. 
However, in the pair plasma case the solution does not 
extend to an infinitely large wavenumber, and there are 
not infinitely many wave modes as we increase u). There 
appears to be a cutoff in frequency, the point where the 
overturn takes place, beyond which there does not exist 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




FIG. 2: Dispersion curves for a relativistic plasma with 7? = 1, 
d)p = 3. Hatted variables are expressed in terms of the non- 
relativistic cyclotron frequency uJc- See text for more detail. 




FIG. 3: Dispersion curves for a relativistic plasma with r? = 5, 
tip = 3. Hatted variables are expressed in terms of the non- 
relativistic cyclotron frequency uic See text for more detail. 

any solution. Secondly, there exists a wave mode below 
the fundamental cyclotron harmonic which was absent 
in the solution provided in [l7j . In fact, we find that 
a solution below the cyclotron fundamental exists for all 
cases, regardless of ry, for = 3 as we show below. This, 
again, is in contrast with the non-relativistic electron-ion 
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FIG. 4: Dispersion curves for a relativistic plasma with 
»7 = 20, (ip = 3. Hatted variables are expressed in terms 
of the non-relativistic cyclotron frequency ljc- See text for 
more detail. 



plasma where there is no wave mode below the first har- 
monic. 

In comparison to the moderately relativistic case, the 
weakly relativistic case is richer in its behavior as well 
as much more structured. The former only has two fre- 
quency modes for a given wavelength, namely a high and 
a low mode, and the latter has many. Moreover, for 
a moderately relativistic pair plasma none of the wave 
modes found in between the two stationary points (dis- 
cussed below) extend to infinitely small wavelengths. In 
fact, there is a limiting wavelength above which these 
modes exist. As reported in [T7| . one can find similar, 
although much less severe, imperfections in the graph- 
ics produced by the contouring algorithm. The culprit 
here is the non linearity of the equation from which the 
solutions are obtained. 



A. Relativistic effects 

As the plasma particles become strongly relativistic, 
with decreasing rj, the dispersion curves undergo drastic 
changes. We now plot all of the dispersion curves shown 
earlier onto a single plot and analyze the progression from 
the weakly relativistic case to the strongly relativistic one 
(see Fig. [5]). The Bernstein waves are strongly absorbed 
near the cyclotron harmonics in the weakly relativistic 
limit where 77 ^ 1. This effect is the strongest at the 
first cyclotron harmonic at which point the phase veloc- 
ity of the wave Vph — > and the wave loses all of its 
energy in heating up the pair plasma. At higher har- 



FIG. 5: Dispersion curves for different values of r] with ujp = 
3. (Solid) The weakly relativistic case with r) = 20, (Dots) 
the mildly relativistic case with r; = 5, and the two strongly 
relativistic cases: (a) (Dash) with rj = 1 (b) (Dot-dash) with 
77 = 0.5 



monies the same phenomenon is repeated, however with 
decreased efficiency. Upon increasing the thermal speed 
of the plasma particles, the cyclotron resonances become 
much less pronounced and start to disappear completely. 
In a hot magnetized pair plasma no resonant interaction 
between the Bernstein wave and the plasma occurs, and 
the solution occupies only a small region of the d) — fcj^ 
space. 



B. Behavior for Small k± 

Unlike the non-relativistic electron Bernstein waves, 
the behavior of the dispersion curves in the limit of van- 
ishing wavenumber is not so straightforward. The hyper- 
geometric function in the integrand of Eq. ((2T|) can also 
be written in the form of an infinite power series Q 

2^3(01, 02; 01,02,03; 2^) = — — TfTT — ^ — (26) 



E 

m=0 



r(ai)r(a2) 
r(ai + TO)r(a2 + to) x"' 
T{bi + TO)r(62 + TO)r(63 + to) to! 



We see that for a; ^ only the m — term has the 
dominant contribution. In that case, we find that the 
hypergeometric function tends to unity. Consequently, 
exx loses its dependence on uj because the only place to 
appears in Eq. pip is inside the hypergeometric func- 
tion. Moreover, this suggests that (jj{k±) is constant in 
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the domain where k± <^ 1 and Eq. (pT|) fails to describe 
the behavior of the dispersion relation for vanishing k±. 
In fact, one has to keep terms in the infinite sum up to 
order m = 2 to obtain any non-trivial solution. However, 
upon doing so one finds that the integral, again, is ex- 
tremely non-trivial and its solution cannot be expressed 
analytically. This is problematic because the dispersion 
solution is no longer a smooth function and may become 
discontinuous for small k±. 

Nevertheless, this problem can be resolved easily. 
Upon cursory inspection of the integrand one finds that it 
is quadratic in kj_ , making it symmetric under the trans- 
formation kj_ —>■ —k±. Furthermore, we demand that 
the dispersion relation be continuous at k± = 0, thus, 
we employ polynomial interpolation to determine the y- 
intercept. In all three figures we use a polynomial of order 
2 or 4, depending on the shape of the curve, to determine 
the stationary points for vanishing wavcnumber. 



this result into the dielectric response function we find 



1_ ^ 



ki [2k^K2iv)JQ 



pe 



X I [cos7r7a) -I- Ho{2k±p)]dp — 1 

Sm TTJUJ J 



(29) 



This integral again has a similar singularity for 70) — n 
for integer n. We can further simplify this equation by 
noting the leading order behavior of the Struve function 
in the limit 6 — > 00, which goes like HQ{b) cx Then, 

in this limit the dominant term in Eq. ([29]) is the cosine 
term. Although this step is not justifiable given the limits 
of integration where the integrand is evaluated for p ^ 1, 
this does not modify the overall behavior of the dielectric 
response function in the large k± limit. Next, we define 



T(?7, t£>) = vro) / dp "fpe '''' cot(7r7cD) 



(30) 



C. Behavior for Large kj_ 

In the weakly relativistic case, one finds wave-particle 
resonances occurring at cyclotron harmonics. These res- 
onances extend to large values of kj_ but not to infinity. 
Intuitively, one may argue, by looking at the rest of the 
dispersion solution, that such a behavior is expected as 
the dispersion curves start and end at fcj^ = 0. The plot 
remains connected over the whole domain, as is evident 
in the strongly relativistic case, and the resonances only 
extend to some maximum value of the wavenumber k^^^. 
This hypothesis can be ascertained by computing the in- 
tegral in Eq. (pij) for fc_L ^ 1, however for large values 
of fcj_ the calculation becomes very computationally ex- 
pensive. 

Ideally, one would like to find the asymptotic behavior 
of the hypergeometric function in Eq. (|2ip to simplify the 
problem. The asymptotic behavior of the hypergeomet- 
ric function can be gleaned by asymptotically expanding 
the Bessel function in the integral representation of the 
hypergeometric function given in Eq. (j20p . The Bessel 
function expansion for large arguments is given as [20j 



^os(.^|.-J)+o(.-='^) (27) 



and solve for exx = 0. After some rearrangement of 
terms, we arrive at a cubic equation which we then solve 
for k± 



2cD ryfcj 



Z(?7,(i) =0 



(31) 



To do the integral in Eq. ((30l) we again employ the same 
contour integration scheme as was done earlier (see Fig. 
[T|). For rj — 20 and C) ~ I, which is the strongest reso- 
nance of all occurring at higher cyclotron harmonics, we 
find a maximum wavenumber k^^^ ~ 72. We find val- 
ues of the same order for resonances at higher harmonics 
as well. This limit shows the maximum value of for 
which a solution exists using the simplified dielectric ten- 
sor (appropriate for large values of A;_l). The actual limits 
to the wavenumber in the realistic dispersion relation are 
typically lower. This exercise demonstrates that the res- 
onances do not extend to infinitely small wavelengths, 
that there is some cutoff at k±_ ~ k^^^, and the dis- 
persion curves remain connected over the whole domain. 
More importantly, one must not forget that this estimate 
is particularly inaccurate for the strongly relativistic case 
where the solution exists for only modest values of k±_ . 



Plugging this back into Eq. (pOI and carrying out the 
integral over the polar angle yields. 



2 

cos(a7r) 



dO cos {b sin 6 
Ho{2b) 



a I cos 

2 4 



/ TT vr 

osinw -I a 

V 2 4 



(28) 



where Hq is the Struve function of order 0. For our spe- 
cific case, b = k±p and a = 71^, and upon substitution of 



D. Stationary Modes 

There are two stationary modes {vg — duj/dkj^ = 0) 
present for vanishing k± in all the cases shown above. We 
plot the evolution of both stationary points as a function 
of 77 in Fig. [6l For Up = 3ujc the upper stationary mode 
remains above the cyclotron fundamental and the lower 
stationary mode remains below it for all 77. Contrastingly, 
in the electron-ion case a solution exists for vanishing k± 
at all cyclotron harmonics, except at the fundamental. It 
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FIG. 6: Stationary modes at vanishing k± as a function of 
non-dimensional reciprocal temperature parameter r/. The 
upper curve corresponds to the upper stationary mode and 
the lower curve to the lower stationary mode. For this plot 
we assume ujp — 3. 

appears that the upper stationary mode turns into the 
hybrid resonance given by 

u;j,=ci,l + l (32) 

in an electron-ion plasma. 

This picture is slightly modified as we lower the plasma 
frequency so that it equals the cyclotron frequency, cDp = 
1, while remaining in the weakly relativistic limit with 
rj — 20 (see Fig. [7]). We still find those two stationary 
points and the strong resonance at the cyclotron fun- 
damental, however, the rest of the plot has disappeared, 
and, interestingly, been replaced by a single closed curve. 
By comparing the present case to the one treated pre- 
viously, with Wp = 3, we find that upon decreasing 
the plasma frequency the lower branch of the dispersion 
curve in the first harmonic band separates from the up- 
per branch. Also the upper branch in the first harmonic 
band connects to the lower branch in the second har- 
monic band. This is the first incidence where a closed 
curve solution, like the ones found by the authors of [l3| , 
has appeared in our analysis. Furthermore, the number 
of stationary modes are now double of what was observed 
previously. 

We also provide a plot of the position of stationary 
points in frequency as a function of the plasma frequency 
for different values of 77 = 1,5,20 (see Fig. [5]). As the 
slope of the curves varies with a change in 77, the up- 
per hybrid frequency, if one was to associate the upper 
stationary mode with it, necessarily depends on the tem- 
perature of the plasma, as shown in Fig. [5] Interestingly, 



FIG. 7: Dispersion curves at two different plasma frequencies 
with 77 = 20. (a) For d)p = 3, and (b) (bold) for cup = 1. 
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FIG. 8: Stationary modes as a function of the plasma fre- 
quency at different values of the reciprocal temperature pa- 
rameter rj. (a) (Solid) For -q — 20, (b) (Dash) for 77 — 10, and 
(c) (Dot) for ?7 = 1. 

in the case of a hot pair plasma (77 < 1) we find that the 
upper stationary mode extends below the cyclotron fun- 
damental for LOp < 1. Thus, an underdense (wp < lUc) 
strongly relativistic pair plasma does not have any Bern- 
stein wave modes above the cyclotron fundamental. An- 
other consequence of this situation is that the two wave 
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modes exist for only small wavenumbers and therefore 
for extremely large wavelengths. We see that for a given 
plasma temperature, the upper stationary mode depends 
linearly on the plasma frequency for cDp > 1. On the other 
hand, the lower mode remains constant. 



IV. DISCUSSION 

In this article, we investigate the behavior of Bernstein 
waves in a uniform, magnetized, relativistic electron- 
positron pair plasma and provide dispersion curves for 
different values of the non-dimensional reciprocal temper- 
ature. The dispersion solutions in all cases are found to 
be remarkably different than the Bernstein modes found 
in the non-relativistic electron-ion case. For a moder- 
ately relativistic pair plasma we find two Bernstein wave 
modes accompanied by two stationary modes for van- 
ishing wavenumber. We do not find closed curve solu- 
tions 17] in all cases but one where the plasma frequency 
equals the cyclotron frequency in the weakly relativistic 
limit. 

As stated earlier, the Bernstein waves in a non- 
relativistic electron-ion plasma propagate undamped in 
the direction orthogonal to the equilibrium magnetic 
field. This might not be true for such waves in a rel- 
ativistic pair plasma. These waves were found to be very 
weakly damped in a weakly relativistic pair plasma [2l| . 
In this article, we only report the real component of Cj as 
the imaginary component is quite non-trivial to calculate 
for the following reason. With Cj complex Eq. ([25]) does 
not hold since the principal value of the integral is no 
longer purely real. As a result, one must calculate the 



residue i?o which itself is a difficult task as the integrand 
in Eq. ([2T|l cannot be easily transformed into the follow- 
ing form where the singularity arises due to a simple pole 
of order m, 



dip) 



{p-poY 



(33) 



and where g{p) is analytic as p ^ po. Nevertheless, we 
do expect very mild damping of the waves (cDi <g; ujr) at 
least for the weakly relativistic case as discussed in |2l| . 
where they find Im(t2;) to be the largest on the upper 
half of the curve that advances towards the cyclotron 
harmonics, and relatively much smaller on the lower half. 
It remains to be seen if damping is at all observed in the 
moderately relativistic case. 

The results presented in this article have important 
implications for all astronomical objects where a magne- 
tized hot pair plasma is present, for example radio pulsars 
[jj , magnetars [22] , pair-instability supernovae [23] . 
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